Genome composition and GC content influence loci distribution in reduced representation genomic studies.

BACKGROUND: Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored. RESULTS: Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected. CONCLUSIONS: Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers.

BACKGROUND: While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts. RESULTS: Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family. CONCLUSIONS: This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Bartonella in bat flies from the Egyptian fruit bat in the Middle East.

In the family of fruit bats, Pteropodidae Gray, 1821, as in the third most diverse group of bats (Chiroptera), the bacterium of the genus Bartonella was detected in several species as well as in a few species of their insect ectoparasites in some tropical and sub-tropical regions of the Old World. The Egyptian fruit bat, Rousettus aegyptiacus (Geoffroy, 1810), is one of the most widespread fruit bats, occurring between South Africa, Senegal, and Pakistan. In this bat species, Candidatus Bartonella rousetti has been detected in three African populations in Nigeria, Kenya, and Zambia. This fruit bat, however, also occurs in the Palaearctic, an area isolating the species geographically and phylogenetically from the Afrotropical part of its distribution range. We screened the blood-sucking bat flies (family Nycteribiidae) from R. aegyptiacus for the presence of the Bartonella bacteria. A rich material of bat fly Eucampsipoda aegyptia (Macquart, 1850), a monoxenous ectoparasite of the Egyptian fruit bats, was collected at 26 localities in seven countries (Egypt, Iran, Jordan, Lebanon, Oman, United Arab Emirates, and Yemen) of the Middle East in 2007-2013. The DNA isolates from the bat flies were subjected to a three-marker (gltA, ssrA, and intergenic spacer region, ITS) multilocus sequence analysis. Based on the amplification of the fragment of ssrA gene by a real-time PCR, 65 E. aegyptia samples from 19 localities in all seven countries were positive for the bacteria. One to five Bartonella-positive individuals of E. aegyptia were collected per one individual of R. aegyptiacus. An analysis of the ITS and gltA genes indicated the presence of an uncultured Bartonella sp., belonging to the Cand. B. rousetti genogroup, identified from populations of the Egyptian fruit bat in Africa. These results support the hypothesis that Bartonella's diversity corresponds to its host's diversity (and phylogenetic structure). Specific lineages of pathogens are present in specific phylogenetic groups of bats.

Identification of 1600 replication origins in S. cerevisiae.

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

Spacer-Programmed Two-Dimensional DNA Origami Assembly.

Two-dimensional (2D) DNA origami assembly represents a powerful approach to the programmable design and construction of advanced 2D materials. Within the context of hybridization-mediated 2D DNA origami assembly, DNA spacers play a pivotal role as essential connectors between sticky-end regions and DNA origami units. Here, we demonstrated that programming the spacer length, which determines the binding radius of DNA origami units, could effectively tune sticky-end hybridization reactions to produce distinct 2D DNA origami arrays. Using DNA-PAINT super-resolution imaging, we unveiled the significant impact of spacer length on the hybridization efficiency of sticky ends for assembling square DNA origami (SDO) units. We also found that the assembly efficiency and pattern diversity of 2D DNA origami assemblies were critically dependent on the spacer length. Remarkably, we realized a near-unity yield of ∼98% for the assembly of SDO trimers and tetramers via this spacer-programmed strategy. At last, we revealed that spacer lengths and thermodynamic fluctuations of SDO are positively correlated, using molecular dynamics simulations. Our study thus paves the way for the precision assembly of DNA nanostructures toward higher complexity.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila.

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

The pan-plastome of Hemerocallis citrina reveals new insights into the genetic diversity and cultivation history of an economically important food plant.

BACKGROUND: Hemerocallis citrina Baroni (Huang hua cai in Chinese) is a perennial herbaceous plant grown for its flower buds that are eaten fresh or dried and is known as the vegetarian three treasures. The nuclear genome of H. citrina has been reported, but the intraspecific variation of the plastome (plastid genome) has not yet been studied. Therefore, the panplastome of this species collected from diverse locations is reported here for the first time. RESULTS: In this study, 65 H. citrina samples were resequenced, de novo assembled, and aligned with the published plastome of H. citrina to resolve the H. citrina panplastome. The sizes of the 65 newly assembled complete plastomes of H. citrina ranged from 156,048 bp to 156,263 bp, and the total GC content ranged from 37.31 to 37.34%. The structure of the complete plastomes showed a typical tetrameric structure, including a large single copy (LSC), a small single copy (SSC), and a pair of inverted repeat regions (IRA and IRB). Many nucleotide variants were identified between plastomes, among which the variants in the intergenic spacer region were the most abundant, with the highest number of variants concentrated in the LSC region. Based on the phylogenetic tree constructed using the ML method, population structure analysis, and principal component analysis (PCA), the panplastome data were subdivided into five genetic clusters. The C5 genetic cluster was mostly represented by samples from Qidong, Hunan Province, while samples from Shanxi and Shaanxi Provinces were classified into the C4 genetic cluster. The greatest genetic diversity was found in the C1 genetic cluster, and the greatest genetic distance between any two clusters was found between the C4 and C5 clusters. CONCLUSION: The resolution of the panplastome and the analysis of the population structure of H. citrina plastomes provide important data for future breeding projects and germplasm preservation.

Mitochondrial genome of Hancornia speciosa gomes: intergenic regions containing retrotransposons and predicted genes.

BACKGROUND: Plant mitochondrial genomes are characterized by high homologous recombination, extensive intergenic spacers, conservation in DNA sequences, and gene content. The Hancornia genus belongs to the Apocynaceae family, with H. speciosa Gomes being the sole species in the genus. It is an siganificant commercial fruit crop; however, only a number of studies have been conducted. In this study, we present the mitochondrial genome of H. speciosa and compare it with other mitochondrial genomes within the Apocynaceae family. METHODS AND RESULTS: A total of 2.8 Gb of Illumina paired-end reads were used to obtain the mitogenome, resulting in 22 contigs that were merged using 6.1 Gb of Illumina mate-pair reads to obtain a circular chromosome. The mitochondrial genome of H. speciosa is circular, containing 63 predicted functional genes, spanning a length of 741,811 bp, with a CG content of 44%. Within the mitogenome, 50 chloroplast DNA sequences, equivalent to 1.72% of the genome, were detected. However, intergenic spaces accounted for 703,139 bp (94.79% of the genome), and 287 genes were predicted, totaling 173,721 bp. CONCLUSION: This suggests the incorporation of nuclear DNA into the mitogenome of H. speciosa and self duplication. Comparative analysis among the mitogenomes in the Apocynaceae family revealed a diversity in the structure mediated by recombination, with similar gene content and large intergenic spaces.

Anopheles arabiensis continues to be the primary vector of Plasmodium falciparum after decades of malaria control in southwestern Ethiopia.

BACKGROUND: Investigating the species distribution and their role in malaria transmission is important as it varies from place to place and is highly needed to design interventions appropriate to the site. The current study aimed to investigate the Anopheles mosquito species distribution and their infection rate in southwestern Ethiopia. METHODS: The study was conducted in 14 malaria-endemic kebeles (the smallest administrative unit), which were situated in eight different malaria-endemic districts and four zones in southwestern Ethiopia. Ten per cent of households in each village were visited to collect adult mosquitoes using Centers for Disease Control and Prevention (CDC) light traps. The larval and pupal collection was done from breeding sites within the villages, and reared to adults. Female mosquitoes were morphologically identified. The head and thorax of adult Anopheles mosquitoes were tested for circumsporozoite proteins (CSPs) using ELISA. At the same time, legs, wings, and abdomen were used to identify sibling species using PCR targeting the rDNA intergenic spacers region for species typing of the Anopheles funestus group and the internal transcribed spacer 2 region genes for Anopheles gambiae complex. RESULTS: A total of 1445 Anopheles mosquitoes comprising eight species were collected. Of 813 An. gambiae complex tested by PCR, 785 (97%) were Anopheles arabiensis, and the remaining 28 (3%) were not amplified. There were 133 An. funestus group captured and tested to identify the species, of which 117 (88%) were positive for Anopheles parensis, and 15 (11%) were not amplified. A single specimen (1%) showed a band with a different base pair length from the known An. funestus group species. Sequencing revealed this was Anopheles sergentii. Among 1399 Anopheles tested for CSPs by ELISA, 5 (0.4%) An. arabiensis were positive for Plasmodium falciparum and a single (0.07%) was positive for Plasmodium vivax. CONCLUSIONS: Anopheles arabiensis continues to play the principal role in malaria transmission despite implementing indoor-based interventions for decades. Sequencing results suggest that An. sergentii was amplified by the An. funestus group primer, producing PCR amplicon size of different length. Therefore, relying solely on amplifying a specific gene of interest in grouping species could be misleading, as different species may share the same gene.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Genomic and cytogenetic analyses reveal satellite repeat signature in allotetraploid okra (Abelmoschus esculentus).

BACKGROUND: Satellite repeats are one of the most rapidly evolving components in eukaryotic genomes and play vital roles in genome regulation, genome evolution, and speciation. As a consequence, the composition, abundance and chromosome distribution of satellite repeats often exhibit variability across various species, genome, and even individual chromosomes. However, we know little about the satellite repeat evolution in allopolyploid genomes. RESULTS: In this study, we investigated the satellite repeat signature in five okra (Abelmoschus esculentus) accessions using genomic and cytogenetic methods. In each of the five accessions, we identified eight satellite repeats, which exhibited a significant level of intraspecific conservation. Through fluorescence in situ hybridization (FISH) experiments, we observed that the satellite repeats generated multiple signals and exhibited variations in copy number across chromosomes. Intriguingly, we found that five satellite repeats were interspersed with centromeric retrotransposons, signifying their involvement in centromeric satellite repeat identity. We confirmed subgenome-biased amplification patterns of these satellite repeats through existing genome assemblies or dual-color FISH, indicating their distinct dynamic evolution in the allotetraploid okra subgenome. Moreover, we observed the presence of multiple chromosomes harboring the 35 S rDNA loci, alongside another chromosomal pair carrying the 5 S rDNA loci in okra using FISH assay. Remarkably, the intensity of 35 S rDNA hybridization signals varied among chromosomes, with the signals predominantly localized within regions of relatively weak DAPI staining, associated with GC-rich heterochromatin regions. Finally, we observed a similar localization pattern between 35 S rDNA and three satellite repeats with high GC content and confirmed their origin in the intergenic spacer region of the 35 S rDNA. CONCLUSIONS: Our findings uncover a unique satellite repeat signature in the allotetraploid okra, contributing to our understanding of the composition, abundance, and chromosomal distribution of satellite repeats in allopolyploid genomes, further enriching our understanding of their evolutionary dynamics in complex allopolyploid genomes.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species.

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee.

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Whole genome sequencing reveals candidate genes involving in PAS resistance in M. Tuberculosis isolated from patients in Thailand.

Para-amino salicylic acid (PAS) was first reported by Lehmann in 1946 and used for tuberculosis treatment. However, due to its adverse effects, it is now used only as a second line anti-tuberculosis drug for treatment of multidrug resistant or extensively drug resistant M. tuberculosis. The structure of PAS is similar to para-amino benzoic acid (pABA), an intermediate metabolite in the folate synthesis pathway. The study has identified mutations in genes in folate pathway and their intergenic regions for their possibilities in responsible for PAS resistance. Genomic DNA from 120 PAS-resistant and 49 PAS-sensitive M. tuberculosis isolated from tuberculosis patients in Thailand were studied by whole genome sequencing. Twelve genes in the folate synthesis pathway were investigated for variants associated with PAS resistance. Fifty-one SNVs were found in nine genes and their intergenic regions (pabC, pabB, folC, ribD, thyX, dfrA, thyA, folK, folP). Functional correlation test confirmed mutations in RibD, ThyX, and ThyA are responsible for PAS resistance. Detection of mutation in thyA, folC, intergenic regions of thyX, ribD, and double deletion of thyA dfrA are proposed for determination of PAS resistant M. tuberculosis.

The Role of Clonal Hematopoiesis of Indeterminant Potential and DNA (Cytosine-5)-Methyltransferase Dysregulation in Pulmonary Arterial Hypertension and Other Cardiovascular Diseases.

DNA methylation is an epigenetic mechanism that regulates gene expression without altering gene sequences in health and disease. DNA methyltransferases (DNMTs) are enzymes responsible for DNA methylation, and their dysregulation is both a pathogenic mechanism of disease and a therapeutic target. DNMTs change gene expression by methylating CpG islands within exonic and intergenic DNA regions, which typically reduces gene transcription. Initially, mutations in the DNMT genes and pathologic DNMT protein expression were found to cause hematologic diseases, like myeloproliferative disease and acute myeloid leukemia, but recently they have been shown to promote cardiovascular diseases, including coronary artery disease and pulmonary hypertension. We reviewed the regulation and functions of DNMTs, with an emphasis on somatic mutations in DNMT3A, a common cause of clonal hematopoiesis of indeterminant potential (CHIP) that may also be involved in the development of pulmonary arterial hypertension (PAH). Accumulation of somatic mutations in DNMT3A and other CHIP genes in hematopoietic cells and cardiovascular tissues creates an inflammatory environment that promotes cardiopulmonary diseases, even in the absence of hematologic disease. This review summarized the current understanding of the roles of DNMTs in maintenance and de novo methylation that contribute to the pathogenesis of cardiovascular diseases, including PAH.

The advantage of intergenic regions as genomic features for machine-learning-based host attribution of Salmonella Typhimurium from the USA.

Salmonella enterica is a taxonomically diverse pathogen with over 2600 serovars associated with a wide variety of animal hosts including humans, other mammals, birds and reptiles. Some serovars are host-specific or host-restricted and cause disease in distinct host species, while others, such as serovar S. Typhimurium (STm), are generalists and have the potential to colonize a wide variety of species. However, even within generalist serovars such as STm it is becoming clear that pathovariants exist that differ in tropism and virulence. Identifying the genetic factors underlying host specificity is complex, but the availability of thousands of genome sequences and advances in machine learning have made it possible to build specific host prediction models to aid outbreak control and predict the human pathogenic potential of isolates from animals and other reservoirs. We have advanced this area by building host-association prediction models trained on a wide range of genomic features and compared them with predictions based on nearest-neighbour phylogeny. SNPs, protein variants (PVs), antimicrobial resistance (AMR) profiles and intergenic regions (IGRs) were extracted from 3883 high-quality STm assemblies collected from humans, swine, bovine and poultry in the USA, and used to construct Random Forest (RF) machine learning models. An additional 244 recent STm assemblies from farm animals were used as a test set for further validation. The models based on PVs and IGRs had the best performance in terms of predicting the host of origin of isolates and outperformed nearest-neighbour phylogenetic host prediction as well as models based on SNPs or AMR data. However, the models did not yield reliable predictions when tested with isolates that were phylogenetically distinct from the training set. The IGR and PV models were often able to differentiate human isolates in clusters where the majority of isolates were from a single animal source. Notably, IGRs were the feature with the best performance across multiple models which may be due to IGRs acting as both a representation of their flanking genes, equivalent to PVs, while also capturing genomic regulatory variation, such as altered promoter regions. The IGR and PV models predict that ~45 % of the human infections with STm in the USA originate from bovine, ~40 % from poultry and ~14.5 % from swine, although sequences of isolates from other sources were not used for training. In summary, the research demonstrates a significant gain in accuracy for models with IGRs and PVs as features compared to SNP-based and core genome phylogeny predictions when applied within the existing population structure. This article contains data hosted by Microreact.

Herbivory induced methylation changes in the Lombardy poplar: A comparison of results obtained by epiGBS and WGBS.

DNA cytosine methylation is an epigenetic mechanism involved in regulation of plant responses to biotic and abiotic stress and its ability to change can vary with the sequence context in which a cytosine appears (CpG, CHG, CHH, where H = Adenine, Thymine, Cytosine). Quantification of DNA methylation in model plant species is frequently addressed by Whole Genome Bisulfite Sequencing (WGBS), which requires a good-quality reference genome. Reduced Representation Bisulfite Sequencing (RRBS) is a cost-effective potential alternative for ecological research with limited genomic resources and large experimental designs. In this study, we provide for the first time a comprehensive comparison between the outputs of RRBS and WGBS to characterize DNA methylation changes in response to a given environmental factor. In particular, we used epiGBS (recently optimized RRBS) and WGBS to assess global and sequence-specific differential methylation after insect and artificial herbivory in clones of Populus nigra cv. 'italica'. We found that, after any of the two herbivory treatments, global methylation percentage increased in CHH, and the shift was detected as statistically significant only by epiGBS. As regards to loci-specific differential methylation induced by herbivory (cytosines in epiGBS and regions in WGBS), both techniques indicated the specificity of the response elicited by insect and artificial herbivory, together with higher frequency of hypo-methylation in CpG and hyper-methylation in CHH. Methylation changes were mainly found in gene bodies and intergenic regions when present at CpG and CHG and in transposable elements and intergenic regions at CHH context. Thus, epiGBS succeeded to characterize global, genome-wide methylation changes in response to herbivory in the Lombardy poplar. Our results support that epiGBS could be particularly useful in large experimental designs aimed to explore epigenetic changes of non-model plant species in response to multiple environmental factors.

Complete chloroplast of four Sanicula taxa (Apiaceae) endemic to China: lights into genome structure, comparative analysis, and phylogenetic relationships.

BACKGROUND: The genus Sanicula comprises ca. 45 taxa, widely distributed from East Asia to North America, which is a taxonomically difficult genus with high medicinal value in Apiaceae. The systematic classification of the genus has been controversial for a long time due to varied characters in key morphological traits. China is one of the most important distributed centers, with ca. 18 species and two varieties. At present, chloroplast genomes are generally considered to be conservative and play an important role in evolutionary relationship study. To investigate the plastome evolution and phylogenetic relationships of Chinese Sanicula, we comprehensively analyzed the structural characteristics of 13 Chinese Sanicula chloroplasts and reconstructed their phylogenetic relationships. RESULTS: In present study, four newly complete chloroplast genome of Sanicula taxa by using Illumina sequencing were reported, with the typical quadripartite structure and 155,396-155,757 bp in size. They encoded 126 genes, including 86 protein-coding genes, 32 tRNA genes and 8 rRNA genes. Genome structure, distributions of SDRs and SSRs, gene content, among Sanicula taxa, were similar. The nineteen intergenic spacers regions, including atpH-atpI, ndhC-trnM, petB-petD, petD-rpoA, petN-psbM, psaJ-rpl33, rbcL-accD, rpoB-trnC, rps16-trnQ, trnE-psbD, trnF-ndhJ, trnH-psbA, trnN-ndhF, trnS-psbZ, trnS-trnR, trnT-trnF, trnV-rps12, ycf3-trnS and ycf4-cemA, and one coding region (ycf1 gene) were the most variable. Results of maximum likelihood analysis based on 79 unique coding genes of 13 Chinese Sanicula samples and two Eryngium (Apiaceae-Saniculoideae) species as outgroup taxa revealed that they divided into four subclades belonged to two clades, and one subclade was consistent with previously traditional Sanicula section of its system. The current classification based on morphology at sect. Sanicla and Sect. Tuberculatae in Chinese Sanicula was not supported by analysis of cp genome phylogeny. CONCLUSIONS: The chloroplast genome structure of Sanicula was similar to other angiosperms and possessed the typical quadripartite structure with the conserved genome arrangement and gene features. However, their size varied owing to expansion/contraction of IR/SC boundaries. The variation of non-coding regions was larger than coding regions of the chloroplast genome. Phylogenetic analysis within these Chinese Sanicula were determined using the 79 unique coding genes. These results could provide important data for systematic, phylogenomic and evolutionary research in the genus for the future studies.

Revisiting gene typing and phylogeny of Trypanosoma cruzi reference strains: Comparison of the relevance of mitochondrial DNA, single-copy nuclear DNA, and the intergenic region of mini-exon gene.

Chagas disease is a widespread neglected disease in Latin America. Trypanosoma cruzi, the causative agent of the disease, is currently subdivided into six DTUs (discrete typing units) named TcI-TcVI, and although no clear association has been found between parasite genetics and different clinical outcomes of the disease or different transmission cycles, genetic characterization of T. cruzi strains remains crucial for integrated epidemiological studies. Numerous markers have been used for this purpose, although without consensus. These include mitochondrial genes, single or multiple-copy nuclear genes, ribosomal RNA genes, and the intergenic region of the repeated mini-exon gene. To increase our knowledge of these gene sequences and their usefulness for strain typing, we sequenced fragments of three mitochondrial genes, nine single-copy nuclear genes, and the repeated intergenic part of the mini-exon gene by Next Generation Sequencing (NGS) on a sample constituted of 16 strains representative of T. cruzi genetic diversity, to which we added the corresponding genetic data of the 38 T. cruzi genomes fully sequenced until 2022. Our results show that single-copy nuclear genes remain the gold standard for characterizing T. cruzi strains; the phylogenetic tree from concatenated genes (3959 bp) confirms the six DTUs previously recognized and provides additional information about the alleles present in the hybrid strains. In the tree built from the three mitochondrial concatenated genes (1274 bp), three main clusters are identified, including one with TcIII, TcIV, TcV, and TcVI DTUs which are not separated. Nevertheless, mitochondrial markers remain necessary for detecting introgression and heteroplasmy. The phylogenetic tree built from the sequence alignment of the repeated mini-exon gene fragment (327 bp) displayed six clusters, but only TcI was associated with a single cluster. The sequences obtained from strains belonging to the other DTUs were scattered into different clusters. Therefore, while the mini-exon marker may bring, for some biological samples, some advantages in terms of sensibility due to its repeated nature, mini-exon sequences must be used with caution and, when possible, avoided for T. cruzi typing and phylogenetic studies.